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We performed numerical simulations of the g-state Potts model to compute the reduced conduc- 
tivity exponent t/v for the critical Coniglio-Klein clusters in two dimensions. At criticality, at least 
for q < 4, the conductivity scales as C(L) ~ L~» , where t and v are, respectively, the conductivity 
and correlation length exponents. For q = 1, 2, 3, and 4, we followed two independent procedures 
to estimate t/v. First, we computed directly the conductivity at criticality and obtained t/v from 
the size dependence. Second, using the relation between conductivity and transport properties, 
we obtained t/v from the diffusion of a random walk on the backbone of the cluster. From both 
methods, we estimated t/v to be 0.986 ± 0.012, 0.877 ± 0.014, 0.785 ± 0.015, and 0.658 ± 0.030, 
for q = 1, 2, 3, and 4, respectively. We also evaluated t/v for non integer values of q and propose 
the following conjecture 40gt/v — 72 + 20g — 3g 2 for the dependence of the reduced conductivity 
exponent on q, in the range < q < 4, where g is the Coulomb gas coupling. 

PACS numbers: 05.50.+q, 64.60.al, 89.75.Da 



I. INTRODUCTION 

The g-state Potts model was initially developed to 
study the ferromagnetically ordered state, but its range 
of applications is much wider and includes, for example, 
problems in materials science [TJ [2] , biology [3J , opinion 
dynamics [1] , image processing [5. , and quantum chromo- 
dynamics [BJ. The general interest for statistical physics 
stems from its rich critical behavior and the fact that 
it generalizes several other models like, e.g., spanning 
trees ((7 = 0), bond percolation (q = 1), and Ising model 
(q = 2) 7J. The critical properties of the order-disorder 
transition in the Potts model depend on the dimensional- 
ity and this transition might even be absent on complex 
geometries [81411] . In particular, in two dimensions, the 
nature of the transition is second order for q < 4 and first 
order for q > 4. 

The theorem of Kasteleyn-Fortuin shows that, for all 
real values of q > 0, the magnetic transition can be de- 
scribed by a purely geometrical model, where the par- 
tition function is a sum over bond percolation configu- 
rations weighted by a factor depending on the number 
of possible states and clusters [T5J. This beautiful re- 
sult has provided the necessary ingredients to develop 
advanced analytic and computational techniques to char- 
acterize the critical properties of the model [131 [T3] . Nu- 
merically, as proposed by Coniglio and Klein, geometrical 
clusters can be obtained starting with magnetic clusters, 
defined as sets of neighboring spins in the same state 
[T5HT7] . and diluting bonds in such a way that, at criti- 
cality, their percolation-like properties are equivalent to 
the magnetic ones [T5H2"U] . 



The geometrical properties of Coniglio-Klein clusters 
have been intensively studied obtaining several numer- 
ical and exact results [5TJ [55] . However, the transport 
properties for q > 1 are, to the best of our knowledge, 
still unexplored. In this work, we study the conductivity 
of these clusters for values of q between zero and four 
in two dimensions. Given the self-similar properties of 
the clusters at criticality for q < 4, we assume that the 
conductivity C between two points scales as, 

C(r)~r-*, (1) 

where r is the distance between the points, t the con- 
ductivity exponent, and v the critical exponent of the 
correlation length, t/v is then the reduced conductivity 
exponent, the focus of our work. 

For bond percolation (q = 1), the conductivity has 
been studied with different methods as, e.g., diffusion 
methods [53J , transfer matrix methods [54T15B] , and star- 
triangle transformations [5TH55] . Some conjectures which 
have been proposed for the conductivity exponent |30[I31] 
were ruled out by numerical data [5BJ. In our study, we 
extend the calculation of the conductivity exponent to 
q = 2, 3, and 4. 

The paper is organized as follows. In Section [H] we de- 
scribe the g-state Potts model and its relation with bond 
percolation. In Section [Hi] we present the two procedures 
to estimate the value of the reduced conductivity expo- 
nent. In Section [TV] we report our results which are then 
discussed in Section [V] The concluding remarks are in 
Section EH 



II. MODEL 
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The Hamiltonian of the model is, 

H = ~ £ J v (<W; - !) > ( 2 ) 

where 5 ai<7j is unity if <7j = Oj and zero otherwise, Jy is 
the coupling strength between spins i and j, and the sum- 
mation is over all pairs of interacting spins. For simplic- 
ity, we assume only interactions among nearest neighbors 
with the same coupling strength, = J. 

A geometrical description of the model which repro- 
duces the scaling behavior of the magnetic transition is 
provided by the Kastelyn-Fortuin theorem 12 . The par- 
tition function associated to the g-state Potts model can 
be written as a sum over all possible bond configurations 

2 = 2^(1 (3) 
M 

where 2Vj is the total number of bonds, and b and N c are, 
respectively, the number of occupied bonds and clusters 
in the configuration v. This is equivalent to a sum over 
configurations of bond percolation in the original graph, 
weighted by a factor of q " , and, consequently, for q = 
1, one recovers the generating function of the random 
bond percolation model. Bonds are established between 
neighboring spins in the same state with probability p = 
1 — e~ K , where K — -^jt, T is the temperature, and ks 
is the Boltzmann constant. Coniglio-Klein clusters are 
then defined as sets of spins (nodes) connected through 
these bonds. This equivalence has been useful in the 
development of efficient Monte-Carlo procedures where 
the critical slowing down is significantly reduced Q21 HI] • 

Let us consider now any two points on the largest clus- 
ter. One can define the backbone as the set of all possi- 
ble paths joining these two points. All the bonds of the 
largest cluster outside the backbone are dangling ends. In 
the backbone, one can distinguish the red from the blue 
bonds: red bonds are singly connected bonds that, if one 
of them is removed, splits the cluster into two, while all 
the others are blue bonds, see Fig. [T] The typical picture 
of the largest cluster at criticality is a set of blue bonds 
forming blobs, i.e., sets of parallel paths, linked together 
by red bonds [32H33- 

Describing the bonds as resistors, one can analyze the 
transport properties of the cluster. In this work, we apply 
a potential difference between two points of the largest 
cluster, separated by a distance r, and we solve the Kirch- 
hoff's law for the network of resistors. We characterize 
the dependence of the conductivity on the distance be- 
tween the points and compute the conductivity exponent 
for different values of q. 

III. SIMULATIONS 

We simulate the g-state Potts model on a square lattice 
with periodic boundary conditions in vertical and hori- 




FIG. 1. (color online) Largest cluster of a random cluster 
configuration for q = 1 and a lattice size of 32. The backbone 
is formed by blobs of bonds belonging to parallel paths (blue 
bonds) linked together by red bonds. The other bonds (yellow 
bonds) are dangling ends. 

zontal directions. The Coniglio-Klein clusters are gen- 
erated using the Swendsen-Wang algorithm [T3]. As de- 
scribed above, critical clusters are identified by establish- 
ing bonds between neighboring spins in the same state, 
with probability p c — 1 — e~ Kc , where K c = k ^ T and T c 
is the critical temperature We focus on the largest 
cluster without periodic boundary conditions. Therefore 
we compute the part of the largest cluster that is in the 
middle of the lattice. In order to do so, we look for a node 
belonging to the largest cluster, starting from the nearest 
point that is down right from the center of the system. 
If it does not belong to the largest cluster, we continue, 
going through the nodes of the lattice by scanning the 
lattice in two directions at the same time: starting from 
the middle point, we go up from left to right and down 
from right to left. When we reach a node that belongs 
to the largest cluster, we clip the fraction of the largest 
cluster, without periodic boundary conditions, contain- 
ing this node. We then identify the lowest and the highest 
node: the lowest node is defined as the first node belong- 
ing to the cluster when one scans through the lattice from 
the bottom to the top (left to right), whereas the highest 
node is the first node belonging to the cluster when scan- 
ning through from the top to the bottom (right to left). 
Since the current vanishes in the dangling ends and only 
flows through the backbone, we restrict our calculation 
to the latter, which can be identified with the burning 
method proposed by Herrmann et al. |36) . Using Kirch- 
hoff 's laws, for every node i in the backbone, we have 

Y^di(V i -V j ) = 0, (4) 

3 

where the sum is over all neighboring spins j, and the 
conductivity C« between nodes i and j, is unity if there 
is a bond between them or zero otherwise. We impose 
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the boundary condition V — L 2 for the highest node and 
zero potential for the lowest one. We then invert the 
(sparse) conductivity matrix, with elements Cy, using 
a sparse matrix solver |37| and obtain the potential for 
every node as well as the global conductivity. Wc com- 
pute the dependence of the conductivity on the distance 
r. For q < 4, we show that the conductivity C as a func- 
tion of the distance r is a power law and we estimate the 
reduced conductivity exponent D a , given by 

C(r)~r- D -. (5) 

The value V — L 2 has been chosen to avoid the sparse 
matrix solver to deal with small potentials, as the back- 
bone mass is scaling with an exponent larger than one 

m 

Simulations were performed on lattices of size 
L = 32, 64, 128, 256, 512, 1024, 2048, for q between one 
and three, and for q — 1 a lattice size of L = 4096 was also 
studied. For q = 4, only system sizes up to 1024 where 
considered. The number of samples generated for each 
lattice size ranges from 10 8 for the smallest system sizes 
till about 10 4 for the largest system sizes. To decorrelate 
from the initial configuration, we reject the first 2L, 3L, 
and 16L sweeps, for q = 2, 3, and 4, respectively. Be- 
tween measurements, we reject eight, 16, and 64 sweeps, 
for q = 2, 3, and 4, respectively. We then use a binning 
procedure and the integrated autocorrelation time from 
Ref. [35] to obtain the error bars. For q = 1, since all 
spins are in the same state, each generated configuration 
is a random bond percolation configuration independent 
from the previous one. 

We also implemented a procedure based on the diffu- 
sion of a random walk on the backbone to estimate the 
reduced conductivity exponent. Such diffusion can be re- 
lated with the transport properties since, on a random 
fractal medium, the walker's mean square displacement 
(R (N)) relates with the number of steps N as 

(R 2 (N))~N^, (6) 

where d w is the random walk dimensionality. The mean 
probability to return to the initial position scales as 

(P Q (N)) ~N-^i ~7V-t, (7) 

where df is the fractal dimension of the backbone and d s 
is the fracton dimension [30] . The reduced conductivity 
exponent can then be computed from the identity [40 , 




The diffusion of a random walk on the backbone can be 
computed exactly following the algorithm proposed in 
Ref. [23]. Once the largest cluster and the lowest and 
highest points have been identified with the method de- 
scribed before, one chooses an initial point in the middle 
of the backbone and restricts the diffusion process to a 



TABLE I. Numerical results for the reduced conductivity ex- 
ponent Da , random walk dimension d w and fracton dimension 
d s . (t/u) IV , is the reduced conductivity exponent obtained 
with Eq. Q . 

1 v _ D ° ^ (£) rw 

1 | 0.987 ±0.013 2.617 ±0.009 1.247 ±0.007 0.985 ± 0.013 

2 1 0.877 ± 0.014 2.587 ±0.012 1.325 ±0.014 0.873 ± 0.022 

3 I 0.785 ± 0.015 2.535 ±0.024 1.396 ±0.025 0.766 ± 0.039 



constant chemical length from this point. To do so, one 
takes the site of the largest cluster which is closer to 
the middle of the line joining the lowest and the highest 
point. The diffusion is solely considered on the subset 
of sites in the largest cluster which are distant from the 
middle point less than a maximum chemical length. This 
distance is defined as the minimum number of steps for 
the walker to go from one site to the other. At each step, 
for every node in the backbone, one computes exactly the 
probability of the walker to be at this node and we mea- 
sure (R 2 (N)) and (Pq(N)). Following this procedure, we 
performed simulations on lattices of lateral size L = 1024, 
averaging results over 10 4 configurations. We limited the 
walker to a maximum chemical distance of 300 steps from 
the origin of the walk. With d w and d s , we deduced the 
value of the reduced conductivity exponent using Eq. pf. 



IV. RESULTS 

We start discussing the cases q = {1,2,3}, followed by 
the special values q = 4 and q = 0. For all cases, we 
compare the estimates for the reduced conductivity ob- 
tained from both the conductivity measurement D a and 
the random walk method (t/v) rw . We use the results 
from both methods to estimate the reduced conductivity 
t/v for different values of q. We also present results on 
the reduced conductivity exponent for non-integer values 
of q, namely, q = 1.5,2.5,3.5. 



A. g = l,2,3 

Figure [2] shows the dependence on the distance r of 
the conductivity obtained with the sparse matrix solver 
method. In the inset, we show the rescaled conductivity, 
C(r)/r~ Da for different values of D a . We obtained the 
values of D a for q — 1, 2, and 3, as summarized in Table 
|IJ where error bars are estimated from the asymptotic 
behavior in the insets of Fig. [2j Figures [3] and [4] show, 
respectively, the data for the mean square displacement 
(R 2 (N)) and the probability of return (P (N)}. Given 
the proposed relations Eqs. Q and ([7]), we obtain d w and 
d s . From Eq. ([8]), we estimate the reduced conductivity 
exponent {t/v) TW . All results are summarized in Table [I] 

The results obtained with the sparse matrix solver 
and with the random walk method are consistent with 
each other. Overlapping the confidence intervals ob- 
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FIG. 2. (color online) Conductivity C as a function of the distance r between the highest and lowest points of the largest 
cluster for (a) 9 = 1, (b) q — 2, and (c) q = 3. The solid lines are guides to the eye of the form C(r) = bcr~ Da . In the inset, 
the rescaled conductivity C(r)/r^ D " is plotted as a function of the distance r for different values of the reduced conductivity 
Ax. 
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FIG. 3. (color online) Mean square displacement of the random walk {R 2 } as a function of the number of steps N for (a) 
9=1, (b) q = 2, and (c) q = 3. The solid lines are guides to the eye of the form {R 2 (N)) = buN^ . In the inset, the rescaled 

2 2 

mean square displacement {R (N))/N d ™ is plotted as a function of the number of steps N for different values of the random 
walk exponent d w . 



tained with both methods, we estimate the value of 
the reduced conductivity to be t/v — 0.986 ± 0.012, 
t/v = 0.877 ± 0.014, and t/v = 0.785 ± 0.015, for q = 1, 
2, and 3, respectively. For q = 1, random bond perco- 
lation is recovered. Our result in this limit is consistent 
with the ones previously reported in the literature using 
different methods 0.9745 ± 0.0015 [26], 0.977 ±0.010 [28], 
0.9826 ± 0.0008(23. 



walk are in agreement. From both methods we estimate 
t/v = 0.658 ± 0.030. 



B. q = 4 

In two dimensions, the nature of the magnetic transi- 
tion in the g-state Potts model crosses over from second 
order, for q < 4, to first order, for q > 4, and, there- 
fore, logarithmic prefactors are typically observed for 
q = 4. For example, several works based on renormaliza- 
tion group theory show that such prefactors are necessary 
for the specific heat, spontaneous magnetization, and sus- 
ceptibility [5TJH2]. We compared, for q = 4, the power- 
law scaling of Eq. ([5| with a power-law including log- 
arithmic prefactors of the form C(r) ~ r~ Da \ \og{r)\ Da . 
We did not observe an improvement of the fit including 
the logarithmic prefactor. Therefore, we used the same 
power-law behavior as considered for q < 4. 

Figure [5] shows the conductivity as a function of r, 
giving D a = 0.658 ± 0.030. From Figs. [6] and [7] we ob- 
tain d w = 2.465 ± 0.020 and d s = 1.467 ± 0.026 and 
using Eq. M (i/^) rw = 0.657 ± 0.037. The result ob- 
tained with the sparse matrix solver and with the random 



C. q = 

As q — > the critical temperature T c diverges and, 
consequently, the critical probability p c vanishes. As a 
result, no criticality is observed at finite temperature for 
q = |43j . However, it is possible to show that in some 
limits, the partition function for q — is equivalent to 
the one of the uniform spanning trees [33], i.e., a sum 
over the ensemble of all possible spanning trees. Here we 
study the conductivity of spanning trees. 

The backbone between any two points on a spanning 
tree is a single fractal path of fractal dimension 5/4 |38j . 
an exponent that can be obtained exactly, for example, 
from the size dependence of the number of red bonds |21l 
[45] . If each bond is a resistor with the same resistance, 
we expect the fractal dimension of the conductivity to 
be the same as the one of the shortest path (and the 
backbone), i.e., D a = 5/4. In Figure [8] we show that this 
result still holds when we consider a spanning tree where 
all the points in the top row (V = L 2 ) are connected with 
all the points in the bottom one (V = 0). 
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FIG. 4. (color online) Mean probability to return to the initial position (Po) as a function of the number of steps N for (a) 
5=1, (b) q = 2, and (c) q = 3. The solid lines are guides to the eye of the form (Po) = bpN a . I n the inset, the rescaled 
return probability (Po(N)} /N 2 is plotted as a function of the number of steps N for different values of the fracton dimension 
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FIG. 5. (color online) Conductivity C as a function of the 
distance r, for q = 4. The line corresponds to the fit with the 
power-law Ansatz, in Eq. ([5]), where — 0.658. In the inset, 
the rescaled conductivity C(r)/r~ Da is plotted as a function 
of the distance r between the highest and the lowest points of 
the largest cluster for different values of D a . 
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FIG. 6. (color online) Mean square displacement of the 
random walk {R 2 (N)) as a function of the number of steps 
N, for q — 4. The solid line is a guide to the eye of the 

form {R 2 (N)} = b R N^ where d w = 2.465. In the inset, the 

2 2 

rescaled mean square displacement (R (N))/N d ™ is plotted 
as a function of the number of steps N for different values of 
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q= 1.5,2.5,3.5 



To obtain intermediary values of the reduced con- 
ductivity exponent for non-integer values of q we have 
computed the conductivity of Coniglio-Klein clusters for 
q = 1.5, 2.5, and 3.5. We generated the bond configu- 
ration using the Chayes-Machta algorithm [?B], which is 
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FIG. 7. (color online) Mean probability to return to the 
initial position (Po) as a function of the number of steps N, 
for q = 4. The solid line is a guide to the eye of the form 
(Po(iV)) = bpN^ where d s = 1.467. In the inset, the 

-d s 

rescaled probability of return (Po(JV)) /N 2 is plotted as a 
function of the number of steps N for different values of d s . 




FIG. 



(color online) Conductivity C as a function of the 



lattice size L, for q — 0. The solid line is a guide to the eye 
of the form C(L) = bcL' ", where D a = 5/4. 

a generalization of the Swendson-Wang algorithm. The 
simulations were done for q — 1.5 and 2.5 for lattice 
sizes between 16 and 1024, and for q = 3.5 between 16 
and 512. Results are averages over 10 5 samples for the 
smallest system sizes to 10 3 for the largest ones. For 
q = 2.5, and 3.5, respectively, we rejected the first 3L 
and 16L sweeps to thermalize and 26 and 358 configura- 
tions between two consecutive measurements [39J. The 
estimates of the reduced conductivity exponent obtained 
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FIG. 9. (color online) Conductivity C as a function of the distance r between the highest and lowest points of the largest 
cluster for (a) q = 1.5, (b) q — 2.5, and (c) q — 3.5. The solid lines are guides to the eye of the form C(r) = bcr~ D " . In 
the inset, the rescaled conductivity C(r)/r~ D ° is plotted as a function of the distance r for different values of the reduced 
conductivity D a . 
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FIG. 10. (color online) Reduced conductivity exponent t/v 
as a function of the Coulomb gas parameter g (data points), 
and conjecture (solid line). 

in Fig. [9] are D a = 0.918 ± 0.013, D a = 0.826 ± 0.015 
and D a = 0.734 ±0.020 for q = 1.5, q = 2.5, and q = 3.5, 
respectively. 



V. DISCUSSION 

In Table [TTJ we summarize the results for the conduc- 
tivity exponents for q < 4. For q > 4, the magnetic 
transition is first order in nature and the backbone (as 
the largest cluster) is compact [35J . Therefore, asymptot- 
ically, no dependence of the conductivity on the distance 
between two points is expected. We have confirmed this 
behavior with simulations for q = 5 (not shown). 

Our estimates of the reduced conductivity are mono- 
tonically decreasing with q. This is in line with the results 
for the dimension of the backbone, which increases with 
q and approaches the spatial dimension for q = 4 |38j . 
As the dimension of the backbone increases, it becomes 
a dense object and the conductivity exponent vanishes, 
being zero for q > 4. 

The conductivity C as a function of the distance r 
shows a bump in the reduced conductivity C(r) /r~ D ° 
around r = 16 for q between one and four. This behavior 
was also previously observed for q = 1 on the square 
lattice, with the star-triangle transformation, and is not 



due to statistical fluctuations [27H29] . 

We have only considered integer values of q in the range 
between zero and four. However, the Potts model can be 
defined for any real value of q [12]. From the obtained 
conductivity exponents we can conjecture a dependence 
on the value of q. The relation between the g-state Potts 
model and the Coulomb gas theory [37] 25] has been very 
useful to develop exact relations for the red bond and 
hull fractal dimensions [3T] H9] . and also to conjecture 
the shortest path fractal dimension [50] , The Coulomb 
gas coupling g is related with q by, 

< ? = 2 + 2cos(^), (9) 

with 2 < g < 4. From the obtained results we propose 
the following conjecture for the dependence of the con- 
ductivity exponent on g, 

£c>-£ + 5-5i* m 

This conjecture fits, within the error bars, the data points 
and displays a continuous behavior in the limit q — > 0, 
where it converges to the value 5/4. It was obtained 
under the assumption that the reduced conductivity ex- 
ponent takes the form £(<?) = ag + b + c/g and that 
it converges towards 5/4 for q — >• 0. The values of the 
constants a, b, and c have been obtained by fitting the 
equation to our results. 

VI. CONCLUSION 

We implemented two independent methods to compute 
the reduced conductivity exponent t/v for different val- 
ues of q in the g-state Potts model. In the first method 
we compute directly the conductivity on Coniglio-Klein 
clusters whereas in the second one we compute exactly 
the diffusion of a random walk on the backbone of these 
clusters. The results obtained with both methods are in 
agreement with each other. From the data of the val- 
ues of the reduced conductivity exponent t/v for integer 
values of q, we propose a conjecture which interpolates 
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TABLE II. Values of the conjectured conductivity exponent 
t/v(g) and the measured reduced conductivity exponent t/u 
for q — 0,1,2,3,4 and the corresponding Coulomb gas pa- 
rameter g. 

g a Ha) j 

2 | | 

1 | 0.975 0.986 ±0.012 
1.5 2.8391 0.921 0.918 ± 0.013 

2 3 0.875 0.877 ±0.014 
2.5 3.1609 0.832 0.826 ± 0.015 

3 f 0.790 0.785 ±0.015 
3.5 3.5399 0.743 0.734 ± 0.020 

4 4 0.650 0.658 ±0.030 
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